Thermodynamic Properties of the One-Dimensional Extended Quantum Compass 

Model in the Presence of a Transverse Field 



R. Jafari 1 - 2 

1 Research Department, Nanosolar System Company (NSS), Zanjan 45158-6591 1, /rarj^] 
Department of Physics, Institute for Advanced Studies in Basic Sciences (IASBS), Zanjan 45137-66731, Iran 

(Dated: January 20, 2013) 

The presence of a quantum critical point can significantly affect the thermodynamic properties 
of a material at finite temperatures. This is reflected, e.g., in the entropy landscape S(T,c) in the 
vicinity of a quantum critical point, yielding particularly strong variations for varying the tuning 
parameter c such as magnetic field. In this work we have studied the thermodynamic properties of 
the quantum compass model in the presence of a transverse field. The specific heat, entropy and 
cooling rate under an adiabatic demagnetization process have been calculated. During an adiabatic 
(de)magnetization process temperature drops in the vicinity of a field-induced zero-temperature 
quantum phase transitions. However close to field-induced quantum phase transitions we observe a 
large magnetocaloric effect. 

PACS numbers: 75.10.Pq; 75.30.Sg; 02.70.-c 



I. INTRODUCTION 

The magnetocaloric effect (MCE), in general, ad- 
dresses the change of temperature of magnetic systems 
under the adiabatic variation of an external magnetic 
field which was discovered in iron by Warburg in 1881 1 . 
The MCE has been widely used for refrigeration due 
to potential room-temperature cooling applications^^. 
For example, adiabatic demagnetization of paramagnetic 
salts was the first method to reach temperatures be- 
low 1KP whereas demagnetization of nuclear spins has 
reached temperatures down to lOOpK^ and is still the 
method of choice for /iK-range cooling^. The cooling rate 
at the adiabatic demagnetization for an ideal para- 
magnet (i.e., a system of non-interacting magnetic mo- 
ments) is equal to T/H, which means linear monotonic 
dependence of temperature on strength of the magnetic 
field. This linear dependence rise out of a direct con- 
sequence of the fact that for any paramagnetic system 
the entropy depends only on the ratio T/H, so for any 
isentrope one gets H/T = const. However, the matter 
could undergo crucial changes in systems with interact- 
ing spins. For instance, in ferromagnets near the Curie 
point one can observe a substantial enhancement of the 
effecP. 

As has been shown in early investigations, quantum an- 
tifcrromagnets are more efficient low-temperature mag- 
netic coolers than ferromagnetic This fact is connected 
with the behavior of the entropy of antiferromagnets. 
The entropy of any antiferromagnet at low temperatures 
displays (at least) one maximum as a function of mag- 
netic field, which usually, according to the third law of 
thermodynamics, falls to zero at T — > 0. 

However, the MCE in quantum spin systems has 
recently attracted scientists attention. From one 
hand, field-induced quantum phase transitions lead to 
universal responses when the applied field is varied 
adiabatically^HS] Q n the other hand, it was observed 



that the MCE is enhanced by geometric frustratiort^M^I, 
promising improved efficiency in low-temperature cool- 
ing applications. More generally, the MCE is particularly 
large in the vicinity of quantum critical points (QCPs). 
The MCE is closely related to the generalized Griineisen 
ratios 
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T {dS/dT) c 



where c is the control parameter governing the quantum 
phase transition which would be external magnetic field 
h for MCE. Using basic thermodynamic relations^, the 
generalized Griineisen ratio Th can be related to the adi- 
abatic cooling rate {dT/dH)^ 
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The magnetic cooling rate is an important quantity for 
the characterization of quantum critical points (QCPs), 
i.e., quantum phase transitions between different mag- 
netic structures under tuning the magnetic field at T — 0. 

Quantum Compass Model (QCM) is simplified model 
which describes the nature of the orbital states in the 
case of a twofold degeneracy^-. Originally, the QCM 
has been used to describe the Mott insulators with or- 
bit degeneracies. It depends on the lattice geometry 
and belongs to the low energy Hamiltonian originated 
from the magnetic interactions in Mott-Hubbard systems 
with the strong spin-orbit coupling^. In QCM the or- 
bital degrees of freedom are represented by pseudospin 
operators and coupled anisotropically in such a way as 
to mimic the competition between orbital orderings in 
different directions. For simplicity, the one-dimensional 
(ID) QCM, is constructed by antiferromagnetic order of 
X and Y ps eudo spin components on odd and even bonds, 
respectivel y 19 ^ 20 !. Moreover, the extended version of the 
ID QCM, is obtained by introducing one more tunable 
parameter, has been studied by Eriksson et o/PS. 
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To the best of our knowledge, thermodynamic proper- 
ties of the extended version of ID QCM in a transverse 
field has not been studied so far. Since the extended 
QCM shows two critical lines (first order and second or- 
der transition lines) , investigation on its thermodynamic 
properties and MCE will be of great importance. How- 
ever, in our recent work, we have shown that extended 
QCM in a transverse field reveals a rich phase diagram 
which includes several critical surfaces depending on ex- 
change couplings 22 . In this paper we exploit the method 
employed in Refs. [5U] and [22] to investigate thermody- 
namic properties of extended QCM in absent and pres- 
ence of a transverse magnetic field. 



where J = ( J\ — J 2 ) - 
andfc = ±^, (j = l,3,- 
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It should be pointed out that though the GS in pe- 
riodic and antiperiodic boundary conditions are slightly 
different in the finite-size system, they are identical in the 
thermodynamic limit and the essential features in finite 
size are also not altered qualitatively. 

Finally, diagonalization is completed by a four- 
dimensional Bogoliubov transformation connecting 
_ k , c?, cP_ u and one thus obtains two different 
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kind of quasiparticles, 



II. HAMILTONIAN AND EXACT SOLUTION 
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Consider the Hamiltonian 



N' 
n=l 
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where J\ and J 2 are the odd bonds exchange couplings, 
Li is the even bond exchange coupling and the parame- 
ters hi and h 2 describe the external magnetic field at the 
even and odd lattice sites respectively which its alternat- 
ing nature is assumed to arise from different magnetic 
moments or different g-factors of the spin. The number 
of sites is N = 2N' and for simplicity we assume periodic 
boundary conditions. The above Hamiltonian (Eq. ([lj) 
can be exactly diagonalized by standard Jordan- Wigner 
transformation 23 24 as defined below, 



2b t b i 



bj = 



which transforms spins into fermionic operators Cj. 

The crucial step is to define independent Majorana 
fermiona^ at site n, = C2 n -i and c£ = c 2n . This 
can be regarded as quasiparticles' spin or as splitting the 
chain into bi-atomic elementary cellars'. 

Substituting for crj, and erj (j = 2n, 2n— 1) in terms 
of Majorana fermions with antiperiodic boundary condi- 
tion (subspace with even number of fermions) followed 
by a Fourier transformation, Hamiltonian Eq. ([!]) (apart 
from additive constant), can be written as 
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where E q k = yj2(a + c) and E p k = ^2(a~ c), c 



b in which 



a = hj + hj + jf + Jj + L\ + 2LxJ 2 cos k, 

b = A^^-h^) 2 + 2J 1 L 2 (J 1 J 2 -h 1 h 2 )cosk + jfL\] 

The ground state (Eq) and the lowest excited state 
(Ee) energies are obtained from Eq.(j3|, 

Eg = -l^iEl + El), E B = -±YX&u-I®> 

k k 

which can be written as a function of a and 6, 



E G = -2j2\/a + Vb, E E = -2j^ \i a - Vb (3) 
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It is clear that the ground state is separated from the 
lowest energy pseudospin excitation by a pseudospin gap 
A = \Ee — Eq\, which vanishes at h\h 2 = J\(J 2 ± L{) 
in the thermodynamic limit. 

It should be mentioned that the exact spectrum and 
the pseudospin gap are the same as that obtained in our 
recent work for h\ = h 2 = /P^l. 

For investigation on the thermodynamic properties of 
the model, we have calculated the free energy per sit. 



/ = -i lirn^ [i InSpeM-PH)] , (4) 

where Sp means a summation over all occupation pat- 
terns. The free energy is easily obtained by using the 
, diagonalized quadratic form of H (Eq.pl)), 
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Entropy (£>), specific heat (C v ) and the normalized 
cooling rate (r c ) are related to free energy via a simple 
thermodynamic relation, 



S = (3 
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Then, from the Eq.|5]) one can easily obtain simple 
analytic expressions for all thermodynamic functions of 
the system. The adiabatic demagnetization curves can 
be found from a direct solution of S(c,T) — const. 



III. EXTENDED QUANTUM COMPASS 
MODEL IN A HOMOGENOUS MAGNETIC 
FIELD (hx = h 2 = h) 

In our recent work, we have investigated the phase di- 
agram of the extended quantum compass model in ho- 
mogenous transverse field by use of the gap analysis 
and universality of derivative of the correlation functions 
(ground statepS. This model is always gapful except 
at the critical surfaces where the energy gap disappears. 
We have obtained the analytic expressions for all critical 
fields which drive quantum phase transitions (QPT) as 
a function of exchange couplings, ho = \J J\(J2 + Li), 
K = \/%Jj2 - Li). 
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FIG. 1: (Color online) Zero-temperature phase diagram of the 
extended quantum compass model in the transverse magnetic 
field for Li = Ji = 1. 



It is useful to recall the zero-temperature phase di- 
agram of the extended compass model in homogenous 
transverse field, see Fig. (uP2i (For simplicity we take 
J\ = L\ = 1). For J2 < 1 and small magnetic field 
h < ho (region (I)) the model is in the spin- flop phase 
(the Need order along the axis (x) which is perpendicu- 
lar to magnetic field is called spin flop). In this region 
tuning the magnetic field forces the system to fall into 
a saturated ferromagnetic (SF) phase above the criti- 
cal field (h c = ho). These states exhibit a gap in the 
excitation spectrum where vanishes at the critical field. 
In the case of J2 > 1 there is antiparallel ordering of 
spin y component on odd bonds under the lower critical 
field h Cl — h n (region (II)) and beyond this critical field 
system goes into a gapped spin- flop phase (region (I)). 
This spin-flop phase exist for h < h C2 = ho- Finally, for 
h > h C2 the ground state is the ferromagnetically polar- 
ized state along the magnetic field (region (III)) which 
exhibits again a gap. These different zero-temperature 
regions, in particular the quantum phase transitions at 
ho and h^ could be reflected by the magnetocaloric prop- 
erties at finite temperature^"^. 

The specific heat of extended quantum compass model 
in transverse field has been plotted in Figs. ^ and Q 
versus temperature (T) and magnetic field (h) for J2 = 2 




FIG. 2: (Color online) The three-dimensional plot of specific 
heat of extended quantum compass model in a transverse field 
versus temperature and magnetic field for J2 = 2. For ex- 
tremely low temperature the maximums of the specific heat 



happen at h ci = 1 and h c 



y/3. 



FIG. 3: (Color online) The specific heat of extended quantum 
compass model in a transverse field versus temperature and 
magnetic field for J2 = 0.8. In this case there is only one 
maximum at h c — VL8 for very low temperature. 




FIG. 4: (Color online) The three-dimensional panorama of 
entropy of extended quantum compass model in a transverse 
field versus temperature and magnetic field for J2 = 2. The 
zero-temperature critical fields are signaled by maximum of 
the entropy at h ci — 1 and h C2 = v3 for low temperature. 



and J2 = 0.8 respectively. As it can be seen from Fig. 
|2]), the specific heat reaches its maximums at QCP for 
extremely low temperatures producing the a small well 
between two critical fields h Cl = = 1 and h c . 2 = ho — 
V3, while there is only one maximum at the critical field 
h c = h = \/L8 in Fig. ^. 



FIG. 5: (Color online) The entropy of extended quantum com- 
pass model in a transverse field versus temperature and mag- 
netic field for J2 = 0.8. In this case the zero-temperature 
critical field (h c = vL8) is specified by maximum of entropy 
for extremely low temperature. 



Figs Q and §5§ show the entropy of the model for 
J 2 = 2 and J 2 = 0.8. They manifest that the maxi- 
mums of the entropy for very low temperature occur at 
the critical fields just like the specific heat. This accumu- 
lation of entropy close QCPs indicates that the system 
is maximally undecided which ground state to chose. As 
the temperature is raised, all the characteristic behaviors 
have been disappeared. 

As it was mentioned earlier, a large MCE also charac- 
terizes a distinctly different class of materials, where the 
low temperature properties are governed by pronounced 
quantum many-body effects. These materials exhibit 
a field-induced QCPs (a zero-temperature phase transi- 
tion), and the MCE has been used to study their quan- 
tum criticality. The adiabatic demagnetization curves of 
extended compass model in the transverse field is pre- 
sented in Figs. ^ and ([7]) for constant initial value of 
magnetic field (hi = 3) and several starting tempera- 
tures (Tj) for J 2 = 2 and J 2 = 0.8, respectively. The 
lowest temperatures of an adiabatic process are reached 
at h\ < h C2 = h and h* 2 > h Cl =h lx . At extremely low 
temperatures the difference between the demagnetization 
fields (h*,h%) and critical fields (h Cl ,h C2 ) becomes very 
small. Clearly, the QCM in transverse field chain cools 
down to lower temperatures near the second critical field 
(ho) than the one close to the first critical field (h n ). 

However, the Griineisen ratio (cooling rate) is notewor- 
thy quantity to specify the second-order transitions since 
it necessarily diverges near QCPs and the divergent be- 
havior obeys the universal scaling law^, 
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FIG. 6: (Color online) Adiabatic demagnetization curves of 
the extended quantum compass model in a transverse field for 
J2 — 2. The lowest temperatures occur at zero-temperature 
critical fields h ci = I and h C2 — 



FIG. 7: (Color online) Lines of constant entropy, i.e., adia- 
batic demagnetization curves of the extended quantum com- 
pass model in a transverse field for J2 = 0.8. The mini- 
mum of isentrop locates at the zero-temperature critical field 

h c = 



r h (T^o,h) = -G hT - i - r , 

n — h c 

where Gh is a universal amplitude and the value Gh = 
— 1 expected for a ^-symmetry in the one dimension^. 
Moreover, the sign of the Griineisen ratio changes as en- 
tropy accumulates near a quantum critical point!^. The 
sign change along with the divergence lead to strong sig- 
natures of the Griineisen parameter near QCPs. The 
cooling rate of the extended quantum compass model is 
plotted versus the magnetic field for J 2 = 2 and J 2 = 0.9 
in Figs. Q and ^ respectively. As it can be seen from 
Fig. ([!]), for J-2 = 2, curves show two quantum phase tran- 
sition at h Cl = hn = 1 and h C2 = ho = \/3 and for 
J2 = 0.9 show just one quantum phase transition at h c — 
ho = vL9. For low temperature, the quantum phase 
transition at h Cl = h n = 1 and h C2 = h = V3 are sig- 
naled by sign changes of the cooling rate Th from negative 



to positive values upon increasing field, see Fig.|8|. How- 
ever in Fig.([9| the cooling rate changes the sign when 
the magnetic field crosses the QCP (h c = ha = vL9). By 
increasing the temperature the magnitude of the peaks 
reduce and all features are washed out which implies 
that thermal fluctuations are already strong enough to 
drive the system to the excited state where no quantum 
phase transition can be seen. In other words, the strong 
enhancement of the magnetocaloric effect arising from 
quantum fluctuations near a /i-induced quantum-critical 
point can be used for finding an efficient and flexible high 
performance magnetic cooling over an extended temper- 
ature range. 
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FIG. 8: (Color online) Magnetic cooling rate (Th) versus mag- 
netic field for different values of temperature for J2 = 2. 
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FIG. 9: (Color online) The variations of cooling rate (F^) with 
magnetic field for different values of temperature for J2 = 0.9. 




FIG. 10: (Color online) The three-dimensional plot of specific 
heat of extended quantum compass model versus temperature 
and the parameter which derives the first-order transition (Ji) 
for J 2 = 2. 




IV. EXTENDED QUANTUM COMPASS 
MODEL (hi =h 2 = 0) 

The complete phase diagram of the extended compass 
model has been reported in Refs. [5T] and They 
have shown that this model is always gapful except at 
the critical line. The multicritical point located on the 
intersection of the first-order (J\/Lx = 0) and second 
order transition (J2/ L\ = 1) lines (For simplicity we take 
Li = l). 

The specific heat of extended quantum compass model 
versus temperature (T) and the parameter deriving the 



FIG. 11: (Color online) The three-dimensional panorama of 
entropy of extended quantum compass model versus temper- 
ature and Ji for J 2 = 2. 



first-order transition (Ji) is plotted in Fig. (10). As it is 
clear, the minimum of the specific heat occurs on the first 
order transition line (Ji = 0) for any arbitrary temper- 
ature. Fig. (11) shows the three-dimensional panorama 



of the entropy versus T and J\ which manifests that the 
maximum of the entropy lies on the first order transition 
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FIG. 12: (Color online) Entropy of the extended quantum 
compass model versus T for J2 = 0.5, J 2 = 0.8, J 2 = 1 and 
J2 = 1.6. Entropy on the gapless line (J2 = 1) shows the 
linear behavior. 



FIG. 13: (Color online) Cooling rate (Fj x ) versus Ji for dif- 
ferent values of temperature, J2 = 2. 



line (Ji =0). It is worthy to mention that the system is 
gapful on this line where the degenerate ground-state sep- 
arates from the excited state. So, the prominent charac- 
teristic of the specific heat and entropy on the first-order 
transition line inherited from the existence of a gap. 

However, the maximum entropy and specific heat falls 
out from the second order transition line (J2 = 1) for 
extremely low temperatures. On this line the system 
is gapless and entropy is linear in T (S oc T) for low 
temperatures while in the gaped cases ( J2 7^ 1) we expect 
activated behavior^ 

S cx exp(-A/T) 

where A is gap in the excitation spectrum. The low 
temperature behavior of the entropy is plotted in Fig. 
(12 1 which verifies the linear and exponential behaviors 



of the entropy on the transition line and gapped phases. 
The asymptotic behavior is S — > for T — > in all cases. 
The cooling rate is plotted versus J\ and J2 in Figs. ( 13 ) 
and (14 1, respectively. 



The cooling rate dependence on J\ is plotted in Fig. 
( 13 1 for different values of temperature. It is shown that 



the first-order transitions are signaled by very sharp and 
pronounced positive and negative peaks at the transition 
point (Ji = 0). Further the only difference between the 
curves, which correspond to different temperatures, is dif- 
ference in the strength of the peaks. Fig. ( 14) shows the 



variations of cooling rate with J2. The QCP pinpointed 
by sign changes of the cooling rate V j 2 from negative to 
positive values upon increasing J 2 . The magnitude of the 
peaks grows rapidly with decreasing the temperature. 




FIG. 14: (Color online) The variations of cooling rate (Tj 2 ) 
with the parameter that derives the second-order transition 
( J2) for different values of temperature, Ji = 1. 



V. SUMMARY AND CONCLUSIONS 

In this paper we have studied the thermodynamic 
properties of the one dimensional extended quantum 
compass model in peresence/absent of transverse field. 
We have presented the specific heat, entropy, adiabatic 
demagnetization curves and Griineisen parameter, which 
is proportional to the magnetocaloric effect, as a func- 
tion of the external magnetic field on the thermody- 
namic limit and at finite temperatures. We have used 
the exact result for the entropy to illustrate that field- 



induced quantum phase transitions give rise to maxima 
of the low-temperature entropy, or equivalently minima 
of the isentropes. This leads to cooling during adiabatic 
(de)magnetization processes where the lowest tempera- 
ture is reached close to the quantum phase transition. 
As a consequence, we found a large positive (negative) 
values of the normalized cooling rate for magnetic fields 
slightly above (below) the critical fields. The general 
features of the entropy should not depend on the specific 
choice of the magnetic field h as control parameter and 



indeed similar behavior is found as a function of the ex- 
change couplings. The low-temperature asymptotics of 
the entropy S is exponentially activated in the gapped 
phases and is linear in T on the second order transition 
line. 
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